Shape resonance for the anisotropic superconducting gaps near a Lifshitz transition: 

the effect of electron hopping between layers 
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The multigap superconductivity modulated by quantum confinement effects in a superlattice of 
quantum wells is presented. Our theoretical BCS approach captures the low-energy physics of a 
shape resonance in the superconducting gaps when the chemical potential is tuned near a Lifshitz 
transition. We focus on the case of weak Cooper-pairing coupling channels and strong pair exchange 
interaction driven by repulsive Coulomb interaction that allows to use the BCS theory in the weak- 
coupling regime neglecting retardation effects like in quantum condensates of ultracold gases. The 
calculated matrix element effects in the pairing interaction are shown to yield a complex physics 
near the particular quantum critical points due to Lifshitz transitions in multigap superconductivity. 
Strong deviations of the ratio 2A/T C from the standard BCS value as a function of the position of 
the chemical potential relative to the Lifshitz transition point measured by the Lifshitz parameter 
are found. The response of the condensate phase to the tuning of the Lifshitz parameter is compared 
with the response of ultracold gases in the BCS-BEC crossover tuned by an external magnetic field. 
The results provide the description of the condensates in this regime where matrix element effects 
play a key role. 

PACS numbers: 74.62.-c,74.70.Ad,74.78.Fk 



I. INTRODUCTION 



For a long time the conventional theoretical mod- 
els for high-temperature superconductivity (HTS) in 
cuprates have been based on a single-component elec- 
tronic system with a single effective band, while com- 
pelling experimental evidence of two bands, resulting 
from different orbitals, crossing the Fermi level was re- 
ported since 1988. 1 Unconventional theories for HTS 
have been proposed focusing on the characteristics of 
multigap superconductivity^— The possible role of a Lif- 
shitz transition-^ [called also electronic topological tran- 
sition (ETT) or quantum phase transitions (QPT) of the 
2.5-order] in HTS was discussed in the frame of the single 
band models In 1993 it was proposed a novel paradigm 
that identifies the quantum mechanism for rising the 
critical temperature in the high-temperature range: the 
shape resonance in superconducting gaps in multiband 
systems where the chemical potential is tuned near a Lif- 
shitz transition in one of the bands£~— 

In the last ten years the experimental research has 
been the driving force for shifting the majority of the 
scientific community toward the novel paradigm that 
HTS emerges in complex systems with multiple electronic 
components in the proximity of a Lifshitz transition. In 
fact two electronic components at the Fermi level have 
been found by many groups in cuprates ^r— showing 
two (or more) strongly correlated bands with tendency 
to phase separation ! 20 ' 21 The two bands in YBaCuO sys- 
tem observed in 1988^ are now confirmed by an increas- 
ing number of experiments^— and recent theoretical 
models consider the new paradigm of a multiband sys- 
tem near a Lifshitz transition! 26 ! 27 The theory of multi- 



band or multigap superconductivity in the clean limit 
was developed since long ago, but it was accepted only 
after the discovery of the two-gap superconductivity at 
high temperature in magnesium diboride WLgB-^ and in 
doped diborides i22r— Later on, low-temperature multi- 
band superconductivity has been found in several ma- 
terials, such as titanates, ruthenates, boro-carbides, se- 
lenides. However, the origin of HTS in MgB2 has been 
assigned by many authors to the strong electron-phonon 
intraband coupling in the a band. On the contrary the 
unique features of multigap superconductivity near a Lif- 
shitz transitio n 29 i 32 i 33 received little attention by the ma- 
jority of the HTS community. High-temperature multi- 
band superconductivity has been discovered in 2008 in Fe 
pnictide multilayers. These materials show a clear weak 
intraband electron-phonon coupling, therefore multigap 
superconductivity at 50K in pnictides has determined in 
these last two years the gradual shift of the scientific opin- 
ion toward the nowadays widely accepted idea that spe- 
cific features of multigap superconductivity may be essen- 
tial for the emergence of HTS<2ir— Recently, the search 
for HTS is focusing on the control of multiband super- 
conductivity in carbon nanotubes— 7 . graphene bilayers^ 
and graphane^ 9 - 

Looking for common features of different high-T c su- 
perconducting materials (e.g., cuprates, diborides and 
pnictides), one could list the following: first, multi- 
gap superconductivity^ - second, multiple electronic 
components^ third, a heterostructure at the atomic 
limit, made of superconducting atomic layers interca- 
lated by spacers made of a different material, forming 
a superlattice^ as shown in Fig. 1; fourth, the tuning of 
the Fermi energy (using any physical, chemical or mate- 
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rial manipulation method, such as doping and superlat- 
tice misfit strain) in order to reach a particular point in 
the electronic band structure.— ~— These are the features 
which were proposed to be essential for the synthesis of 
new high-T c superconductors, with ^.eventually reach- 
ing room temperature, in Refs. la-llG. The particular 
condition yielding high T c was thereby identified as the 
tuning of the Fermi energy at a shape resonance for the 
superconducting gaps near a Lifshitz transition. In this 
work we first discuss the fundamental points of shape 
resonances in multiband superconductors made of a first 
two-dimensional (2D) cylindrical Fermi Surface and a 
second Fermi surface, tuning the chemical potention near 
the Lifshitz transition. We present the results for the a 
system of weak intraband coupling so the system is well 
described by the BCS weak-coupling limit and we con- 
sider a case where the interband pairing is the dominant 
interaction so that it could describe the case of pnictides 
where the interband pairing is the dominant interaction. 
In this regime we investigate the effect of variable electron 
hopping between the layers that is the characteristic fea- 
ture of superlattices of superconducting layers at atomic 
limit. 




FIG. 1. (Color online). Pictorial view of artificial heterostruc- 
tures realized as superlattices of stacked planes, with partic- 
ulars of the superconducting layer A and of the spacer layers, 
in four different families of compounds: (a) cuprates, (b) 
magnesium diborides, (c) iron pnictides, and (d) graphene 
or graphane layers. 



II. SHAPE RESONANCES FOR 
SUPERCONDUCTING GAPS 

The concept of shape resonance for superconducting 
gap parameters was first introduced by Blatt and Thomp- 
son for a single 2D membrane^ 6 - and it has been de- 
veloped for a single nanowire.— The shape resonance 
for the superconducting gaps in superlattices of quan- 
tum wells, (quantum layers, quantum wires and quantum 
dots)^^ has been proposed to yield three-dimensional 
(3D) superconductors. Shape resonances in a superlat- 
tice of stripes in cuprates ^r— and pnictides^ 6 ., of layers 
in diborides ) 33 ' 51 " — and in a superlattice of nanowires 
in carbon nanotubes, have been discussed.— The novel 
scenario for a multiband system with a dominant role of 
interband pairing and negligible intraband coupling has 
been proposed before the discovery of pnictides. 

The definition of shape resonance in the theory of 
superconductivity^ 6 - is borrowed from nuclear physics. 
Ettore Majorana first in 1929-1930 developed a theo- 
retical model to describe the artificial disintegration of 
nuclei by bombardment with a-particles^ By follow- 
ing the quantum dynamics of a state resulting from 
the superposition of a discrete state with a continuum 
one, whose interaction is described by a given poten- 
tial term, Majorana has made the first application of 
the concept of the idea of quasi-stationary states inter- 
acting with the continuum j 55 ' 57 The theory of shape- 
elastic scattering cross section for nucleon-nucleus res- 
onances, where the total energy of incoming or outgo- 
ing particles and the quasi stationary states of the com- 
pound nucleus are nearly degenerate, was developed in 
the fifties by Blatt, Feshbach and de Shalit^§2, The 



shape resonance in nuclear scattering from a potential 
exhibits characteristic peaks, as a function of energy, 
for values of energy such that an integer number of 
wavelengths fits within the nuclear potential well. The 
shape resonances have been measured in electron-atom 
scattering 6 ^ and electron-molecule scattering ] 62 ' 63 The 
name Feshbach resonance was coined to indicate a par- 
ticular shape resonance, where the quasi-bound state is 
near zero energy^ The Feshbach resonance in atomic as- 
sociation and dissociation processes occurs in ultra-dilute 
and ultra-cold atomic gases where the energy of a di- 
atomic molecule (tuned by an external magnetic field) 
is degenerate with the chemical potential of the atomic 
gas^ In this condition the exchange interaction give 
the Feshbach resonance. The discovery that in boson 
and fermion ultracold gases the energy of the diatomic 
molecule can be tuned (by means of a magnetic field) 
above and below the continuum threshold has driven to 
the realization of the Bose-Einstein condensation (BEC) 
in boson gases and of the BCS-BEC crossover in fermion 
gases driven by the Feshbach resonance! 6 -^— Therefore 
the shape resonance and Feshbach resonance belong to 
the class of Majorana-Fano-Feshbach resonances due to 
configuration interaction effects between open and closed 
scattering channels described by the Feshbach-Fano par- 
titioning method) 6 ^— In a multiband superconducting 
system, when the chemical potential is tuned at the Lif- 
shitz transition in one of the bands, the pairs of elec- 
trons at the Lifshitz transitions arc in a quantum crit- 
ical point where the group velocity of electrons goes to 
zero and they form bosonic pairs that can be considered 
the equivalent of the Majorana quasi-stationary states or 
the diatomic molecule in the magnetic field in ultra cold 
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gases. 

In superlattices, the shape resonance^— in the 
superconducting gaps appears when the Fermi energy 
of electrons in a band is tuned around the Lifshitz 
transitions^ in a different band. There are several types 
of Lifshitz transitions: 

type I) when the Fermi energy crosses the band edge, 
at E e d ge , with the appearance or disappearance of a new 
Fermi surface (FS); 

type II) when the Fermi energy crosses the electronic 
topological transition ETT^,o-2D where one FS changes 
from 3D to 2D (i.e., from "spherical" to "cylindrical" ) 
or vice versa, with the opening or closing of a neck in a 
tubular FS; 

type III) at the singular point where the FS changes 
from 2D to one-dimensional (ID) or vice versa, with the 
change of FS topology from a closed FS circle to dis- 
connected Fermi arcs^ Once the heterostructure of ma- 
terials selected as the building blocks for the supercon- 
ducting units and for the spacer layers has been synthe- 
sized with its structural parameters, the tuning of the 
Fermi level at an ETT can be controlled by means of: i) 
the "charge transfer" between the superconducting layers 
and the "spacers"; ii) the superlattice misfit strain^— 
between the superconducting and spacer layers; iii) the 
thickness of the spacers; iv) the ordering of dopants in 
the spacer; v) the superstructures in the superconducting 
layers forming stripes; vi) pressure. 



III. BREAKDOWN OF THE STANDARD BCS 
APPROXIMATIONS 

The heterostructures at the atomic limit of supercon- 
ducting units, where the chemical potential is tuned at a 
type I Lifshitz transition (appearing or disappering of a 
new FS spot) and a the type II Lifshitz transition (open- 
ing a neck in a FS), where the superconducting gaps show 
a shape resonance, have been called super stripes 
From the theoretical point of view, the shape resonance 
implies the breakdown of two main approximations of the 
standard BCS theory for superconductivity: 1) the infi- 
nite distance of the Fermi energy from the band edges, 
which is reasonable whenever the energy cutoff for the 
pairing interaction, luq, is much smaller than the Fermi 
energy Ep (as measured from the nearest singular point 
in the electron spectrum): in this case the electron den- 
sity of states (DOS) can be taken as nearly constant 
within an energy window of the order of luq; 2) the single 
band approximation, which is altogether reasonable in 
the dirty limit, when impurity scattering mixes the vari- 
ous electron components yielding an effective single-band 
system. 

The situation in which the first approximation breaks 
down (e.g., when the Fermi energy falls near an 
electronic topological transition) has been called the 
van Hove scenario^ - — or the Pomeranchuk instability 
scenario^Sr— belonging to the class of quantum phase 
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FIG. 2. Pictorial view of the evolution of the FS of the two- 
band electronic system near a band edge crossing two Lif- 
shitz transitions. The chemical potential is tuned so that Ef 
crosses the band edge E e d g e of the second band (the type I 
Lifshitz transition) and the 3D-2D ETT at Ezd-id (the type 
II Lifshitz transition). Here, £ is the transversal band dis- 
persion. On the left side a first large 2D FS coexists with a 
second tubular 2D FS; with decreasing Ef and crossing the 
energy Ezd-2d where the second FS undergoes a 3D-2D ETT 
changing its topology, the tubular 2D FS becomes a closed 3D 
FS as shown in the central panel; only the first large 2D FS re- 
mains when the chemical potential is moved below the band 
edge energy E e d g e of the second band as shown in the left 
panel. 



transition scenarios investigated for pairing in single band 
superconductivity^ - — While these scenarios have been 
proposed for cuprates assuming a single-band approxima- 
tion, recently there is increasing evidence for a multiband 
superconductivity also in cuprates^^ ~ 20 ' 22 ~ 28 ' 41 i 89 ' 90 

As far as the second approximation is concerned, as 
we have recalled above, the very concept of multiband 
superconductivity is meaningful only in the clean limit. 
Indeed, in nearly all metals and alloys several bands 
cross the Fermi level, but impurity scattering leads to 
a mixing of the electron states. Therefore, the single- 
band description, which is adopted, e.g., in the standard 
BCS theory, is a reasonable approximation in the dirty 
limit. However, several experiments in doped diborides 
have shown that the clean limit is robust in heterostruc- 
tures at the atomic limit where both the different parity 
of the bands and their different spatial location forbid 
the mixing of the different electron states at the Fermi 
surfaces.— Nowadays it is commonly accepted that the 
dirty-limit approximation breaks down in magnesium di- 
boride, in pnicitdes, in superlattices of carbon nanotubes, 
in graphene and graphane layers, and multigap supercon- 
ductivity may emerge also in the presence of lattice disor- 
der due to dopants and misfit strain of the superlattice, 
when these are not apt to mix the different electronic 
components. 

Another crucial feature of multiband superconduc- 
tivity is that, beyond the standard intraband attrac- 
tive interactions that promote pairing within each band, 
exchange-like interband interactions, that scatter Coop- 
ers pairs from one band to another, become relevant. In 
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the hypothetical absence of interband interactions, each 
band would be characterized by its own superconduct- 
ing critical temperature. Arbitrarily weak interband in- 
teractions lead to a single critical temperature, like in 
the proximity effect. It is important to notice that in- 
terband interactions are generically repulsive and in this 
case interband pairing leads to condensate wavefunction 
with opposite signs (the so called s± pairing) i n i 34 
In the shape resonance scenario for multiband supercon- 
ductivity, the control of the ratio between the intensity 
of exchange-like interband pairing and intraband Cooper 
pairing, by material design techniques, is crucial. 



IV. TWO BAND SUPERCONDUCTIVITY AT A 
BAND EDGE 

In this work we study the shape resonance in the case 
when the Fermi energy is tuned at a band edge energy 
E e d ge in one of the bands in a multiband superconductor. 
In particular we focus on superlattices of metallic layers, 
with finite hopping probability between the metallic lay- 
ers. The universal feature of band edges in a superlattice 
of metallic layers is that the small FS appearing or dis- 
appearing as the Fermi energy crosses the band edge has 
a 3D topology, due to finite electron hopping between 
the layers separated by a finite potential barrier. We 
consider the interesting case where the electron hopping 
between layers is so small that the transversal band dis- 
persion £ due to hopping between layers is much smaller 
than the energy separation between the subbands. In 
these superlattices the FS appearing or disappearing as 
the Fermi energy crosses the band edge has a 3D topol- 
ogy in the range < Ep — E e dge < £• When the Fermi 
energy crosses the van Hove singularity, an ETT3D-2D 
of the Lifshitz type called "opening a neck" takes place, 
with the FS changing topology, from closed to tubular, 

at Ep — E e dge = £• 

It has been show n n i 36 that the optimum amplifica- 
tion of the critical temperature occurs when £ is of the 
order of the energy cut-off of the interaction, ujq. In 
this case, in the narrow energy range 3wo where the 
chemical potential crosses both the band edge and the 
ETTzd~2D 1 the condensate in the newly appearing FS 
undergoes a crossover from a mixed Bose-Fermi regime 
to the Fermi BCS regime. For example, in the regime 
— 1 < (Ep — Eedge) /^o < 0, the electron states asso- 
ciated with the newly appearing FS are unoccupied in 
the normal state, and a BEC-like condensate of pairs 
with bosonic character is formed below the critical tem- 
perature. In the range < (Ep — E e dge)/^o < 1, ah 
few electrons in the newly appearing 3D FS in the nor- 
mal state, condense below T c , with the breakdown of the 
standard approximation (Ep — E e dge)/^a ^ 1. In the 
range < (Ep — E e d ge )/uio < 1, in the case under dis- 
cussion £ = ljo, the new FS has a 3D topology. The 
ETT3D-2D plays a key role in the configuration interac- 
tion between pairing channels in different bands in the 



range < (Ep - E ed ge)/^o < 2. 

We are interested in the evolution of superconducting 
gaps in the two bands and of the critical temperature in 
the narrow energy range where the Feshbach-like shape 
resonance in superconducting gaps takes place. So, we 
consider a first large 2D FS in the standard BCS ap- 
proximation (i.e., where the standard BCS approxima- 
tion is valid, large Ep/ujQ ratio, i.e., Ep far from band 
edges) that coexists with a second appearing small FS, 
where the standard approximation breaks down (small 
Ep/ojQ ratio, with Ep close to both the band edge and 
the ETT3D-2D): see Fig. 2. 

This simple model is apt to capture the physics emerg- 
ing from experiments in doped diborides^ 3 - cuprateS ) 26 i 27 
and pnictidesi2i~— Moreover it represents the optimal 
case to obtain the enhancement of critical temperature 
via the shape resonance mechanism^ It describes the 
multiband pairing at the bottom or top of 2D bands in 
a generic superlattice of layers. Here, the effective po- 
tential barrier and the thickness of the spacer layers are 
selected to adjust the electron hopping between super- 
conducting layers and the transversal energy dispersion 
£ of the order of the energy cut-off ujq of the interaction. 

The shape resonance is determined by the relative 
strength of the intraband coupling constants and the in- 
terband coupling constant determined by the type of ma- 
terial forming the superconducting layers and the type 
of material forming the spacers. We consider here the 
case, typical of pnictides, where the intraband attractive 
coupling strength is weak in both the first and second 
subband. We have obtained the evolution of the su- 
perconducting gaps from below the edge, in the mixed 
Fermi-Bose regime, to above the ETT in the Fermi-like 
regime. The direct evidence for the quantum interference 
effect between pairing channels is provided by minima 
in the gap parameter for electrons in the large FS. We 
report the evolution of this scenario with increasing ex- 
change interband pairing. We also investigate the effect 
of the variation of the transversal dispersion £ that can 
be changed by changing the separation space between 
layers. These results can explain the differences among 
1111 pnictides, with changing the rare earth ionic ra- 
dius, and the difference between 1111 and 122 pnictides, 
where the spacer layers changes in the superlattice, while 
the superconducting layers are kept fixed. 

Finally, we show that the plots of the BCS gap ratios 
versus the critical temperature that display significant 
changes in the different superconducting regimes of the 
Feshbach-like shape resonance. These plots can be eas- 
ily compared with experimental data and provide a very 
good experimental test for the theory of the Feshbach-like 
shape resonance. 

There is a clear analogy between the Bose-Fermi 
crossover case studied here and the BEC-BCS crossover 
in ultracold Fermi gases.— The theory presented here fol- 
lows Blatt's approach^, and it is appropriate for dealing 
with the case of a chemical potential near a band edge 
since the equation for the gaps is solved together with 
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FIG. 3. (Color online). The shape resonance occurs by tun- 
ing the Fermi energy near the band edge E e d ge of a second 
electron-like band (upper panel) or a second hole-like band ( 
lower panel) . The second band coexists with a large first band 
(with the Fermi energy far from its edges) that has a free- 
electron like dispersion. The second band has a free-electron 
dispersion in the direction of the electron wavevector parallel 
to the plane (fc//) while in the perpendicular direction (k z ) 
of the superconducting layers it is determined by a periodic 
potential barrier determined by the superlattice. This second 
band has a 3D character in the energy range between the band 
edge, and the van Hove singularity energy Eett 3D _ 2D ■ The 
band dispersion is £=4t z , where t z is the electron hopping 
integral between the layers. 



the equation for the chemical potential in the supcrfluid 
phase. In fact like for the BCS-BEC crossover, as shown 
by Legget1j22 the BCS wave function corresponding to an 
ensemble of overlapping Cooper pairs at weak coupling 
considering a contact interaction (BCS regime) evolves to 
non-overlapping pairs with bosonic character as the den- 
sity decreases at the band edge (BEC regime). The BCS 
theory remains valid in this limit if the BCS equation for 
the gap is coupled to the equation that fixes the fermion 
density so that the chemical potential fi results strongly 
renormalized below the critical temperature with respect 
to the Fermi energy Ep of the non interacting system. 

In ultracold gases the energy of the bound state of the 
diatomic molecule above or below the continuum is tuned 
by a magnetic field. In the present case of superlattices 
the chemical potential can be tuned near the ETT of a 
narrow band, for example by superlative misfit strain. 



V. THE TWO BAND MODEL IN A 
SUPERLATTICE OF METALLIC LAYERS 

As pointed out above, the standard multiband BCS 
theory in the weak-coupling limit (or the Eliashberg the- 
ory in the strong-coupling limit )2r— assumes a DOS which 
is (at least approximately) constant within a window of 
thickness 2luq around the Fermi energy. This assumption 
is certainly not valid if the chemical potential is located 
near a band edge or close to a van Hove singularity of a 
superlattice of metallic layers. 

Near the band edges in the l-th band, where Vk-Ec = 0, 
the energy of an electron can be approximated by a free 
electron dispersion law -E^k = El + (k 2 /2m^), where mi 
is the electron effective mass at the band edge. Therefore, 
one can use the free electron approximation in a narrow 
region around the band edge, in the case when the band 
width is much larger than 2wo, the energy range of in- 
terest for pairing. In the case of interest here, however, 
while this approximation is valid in the x, y plane of the 
superconducting layers, it is certainly not valid in the z 
direction, where the dispersion is on the order of ljq. In 
this case, an anisotropic band with weak dispersion in the 
z direction and larger dispersion in the x, y plane above 
the band edge of the second band should be adopted, 

E™=E 2 +E(k z ) + %±^ (1) 

where E(k z ) is the actual energy dispersion in the peri- 
odic potential of the superlattice and mi i is the effective 
mass in the x, y plane. This situation is obtained within a 
model where a free electron gas is confined in a potential 
which is periodic in the z direction, 

n— +oo 

W(z) = W b (z-nd), (2) 

n— — oo 

where W b (z) = -V b for \z\ < L/2 and W b (z) = for 
L/2 < \z\ < d/2, L is the width of the confining well and 
d is the periodicity of the superlattice in the z direction. 
This periodic potential mimics the phenomenology, e.g., 
of the pnictides, diborides, and stacks of graphene layers 
made of a superlattice of stacked planes, as shown in Fig. 
1. The confining potential generates a band structure or- 
ganized in subbands. Each subband has a dispersion in 
the transversal direction, as shown in Fig. 2, that gives a 
Fermi surface with 3D topology, with closed isoenergetic 
surface, near the lower band edge E e d ge — Eg^, and a 
2D character, with isoenergetic surfaces open in the z di- 
rection, above some energy threshold Eett zd -id — Et,T- 
The present model is apt to describe quantum interfer- 
ence phenomena between different scattering channels in 
a large and a small FS which are the object of this work, 
when the chemical potential /i of the system is tuned 
near the bottom of a £-th subband, within a window of 
width 4ujq. In the energy window of width 2wo the first 
band, with a large 2D tubular FS, has a constant DOS 
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N%. The second FS appears as the Fermi energy crosses 
the level E edge = E e ^ L , and changes from closed 3D to 
tubular 2D topology as the Fermi energy crosses the level 
Eett 3D - 2 d — Ei,T as it is shown in Fig. 2. The model 
for electron-like FS can be easily extended to hole-like FS 
as shown in Fig. 3. The DOS of the second FS N2 near 
the edge has the typical 3D behavior. 

In order to solve the BCS equations in the shape res- 
onance scenario it is necessary to determine the electron 
wavefunction. Using the model of a free electron gas 
confined in a superlattice, the wavefunction of electronic 
states can be obtained by solving the Schroedinger equa- 
tion, so that we can calculate the anisotropic k- depen- 
dent gaps in the wavevector space and the interference 
between the intra-band pairing channels and the pair ex- 
change determined by the interband pairing terms. The 
solution of the Schroedinger equation for the ID periodic 
potential of the superlattice of layers allows us to cal- 
culate the gaps also in correspondence of the 3D to 2D 
electronic transition ETT3D-2D at the van Hove singu- 
larity, beyond the standard BCS approximations . 46 i 50 

Going from a single slab^ to a 3D superlattice^, a 3D 
condensate is formed, reducing the effect of fluctuations 
of the superconducting order parameter, which suppress 
the mean- field superconducting T c . The crossover from 
2D to 3D can be described in our model by changing 
from an infinite potential barrier between the planes of 
the superlattice in the z direction, provided by spacer lay- 
ers, to a finite barrier This yields a finite hopping 
term that broadens the sharp discontinuity of the DOS 
of a pure 2D band. This broadening increases the width 
of the shape resonance and, at the same time, yields a 
3D condensate. Moreover, it is possible to design artifi- 
cial superlattice heterostructures at optimum shape res- 
onance condition, i.e., where the value of the potential 
barrier V& and its width are such that the subband dis- 
persion in z direction £ is of the order of energy cut-off 
cjo of the interaction. Therefore, it is possible to tune 
the system to the optimum shape resonance condition, 
achieving an enhancement of the critical temperature. 

The intrinsic k-dependence of the pairing interaction 
k' m t ne superlattice with wave vector k z induces a 
structure in the k-dependent interband coupling interac- 
tion for the electrons that determines the quantum inter- 
ference between electron pairs wavefunction in different 
subbands of the superlattice^ We calculate the term 
V k ' k / following our previous work for a superlattice of 
wires^ - determining the matrix elements of local inter- 
action potential and then introducing a cutoff with two 
9 functions, 
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where k = k z (k' = k z ) is the superlattice wavevector, in 
the z direction, perpendicular to the planes, of the initial 
(final) state in the pairing process, and 
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where Nq(Ef) is the DOS at Ep for a free electron 3D 
system, V^d is the volume of the system, 
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and the V^,k(z) are the eigenfunctions in the superlattice 
of quantum wells, normalized so that J d dz\il)e t k(z)\ 2 = 1. 
The use of single cut off in two band superconductivity 
has been justified in detail in refj^SS 

The dimensionless factor ci y = {— l) Se - e 'c® £l assumes 
positive values for 1 = 1' (intraband Cooper pairing) and 
negative values for £ ^ £' (repulsive exchange-like inter- 
band pairing, with ci y i> = cg'j) and measures the relative 
intensity of intraband and interband pairing strength. In 
fact, it multiplicates the k-dependent integral and there- 
fore permits to simulate the behavior of different super- 
conductive multilayer compounds controlling the ratio 
between intensities of intraband and interband pairings. 

In order to determine the gaps and chemical potential 
self-consistenly and to calculate the superconductive T c 
we use iterative solving methods for the coupled BCS-like 
equations 



M *t 2J(Et>p - p) 2 + Af, ik , 



-y 

72 



d 2 



Eim — p 



{E tM - (j,) 2 + A; 



.k 



(6) 



(7) 



starting with an initial gap parameter equal to a constant 
and an initial chemical potential equal to the Fermi level 
in the normal state and considering the convergence oc- 
curred for relative variation of the gap and charge density 
p less than 1CP 6 . Here M is the total number of wavevec- 
tors k' and p is the electron density. The superconducting 
critical temperature T c is calculated by iterativly solving 
the linearized equation 



A 1 V V"' tanh (jjg 



.k' 



(8) 



until the vanishing solution is reached with increasing 
temperature. Here, the gaps depend on the superlattice 
wavevector k as well as on the subband index. Hence, T c 
and the gap at a given point of k-space become implicit 
function of all the different values of the gaps in the entire 
k-space. 

In the standard BCS theory, where the Fermi energy 
is far from the band edges, the relative variation of the 
chemical potential going from the normal to the super- 
conducting state is expected to be negligible. This is 
not true when the chemical potential is tuned near the 
band edge of the second band. In fact, our calculation 
yields a significant variation of the chemical potential in 
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FIG. 4. (Color online). The case of weak coupling in the sec- 
ond band (the so called pnictide case) where the intraband 
coupling parameters ratio is fixed at C2,2/ci,i = 0.45. The 
critical temperature T c (panel a), the ratio 2Ai/T c (panel 
b) and the ratio —2/S.2/T c (panel c) are plotted as function 
of the reduced Lifshitz parameter (p, — E2.l)/^o- The dif- 
ferent curves in each panel represent the cases of different 
interband pairing strength: the case for C12/C11 = —2.73 
(violet curve with open triangles); the case for C12/C11 = 
— 1.59 (light blue curve with open triangles); the case for 
C12/C11 = —1.04 (green curve with open diamonds); the case 
for C12/C11 = —0.91 (dark blue curve with open squares) and 
shows the case for C12/C11 = —0.68 (red curve with open cir- 
cles). For the largest case of interband pairing strength ci,2 
considered here, the critical temperature reaches 300K in the 
range 1 < (fi — E2,l)/ojo < 2. The shape resonance antireso- 
nance appears as a minimum in the critical temperature for 
(fi — E2,l)/uq=-1 for higher interband repulsive interaction 
and moves towards (/j, — E2,l)/uq = for the smaller inter- 
band repulsive interaction. The ratio 2Ai/T c for the first 
band shows maxima and minima like the critical tempera- 
ture. On the contrary the ratio 2A2/T C shows a minimum 
at (n — E2,l)/^o = and it has a maximum at the value 
of maximum T c for weak interband repulsive coupling, on the 
contrary it exhibits a maximum in the Bose-like regime(yellow 
region) (0 > fi — E2,l)/ojo > — 1 and a minimum value in cor- 
respondence of the maximum T c for very strong interband 
repulsive pairing. 



the superconducting phase, as a function of the charge 
density. A relative variation of the chemical potential 
going from the normal to the superconducting phase, as 
large as 10 -3 , is obtained near the band edge and at the 
ETT30-2D, within a range 4wo- The variation starts to 
be large, as compared with the standard BCS result, in 
proximity of the Bose-Fermi crossover regime below the 
band edge up to well beyond the ETT^jj^d- 

The Fcshbach-like shape resonance regime occurs prop- 
erly in correspondence of this large variation of the chem- 
ical potential between the normal and superconducting 
phases. Our theoretical approach provides a direct mea- 
sure of the gap, at a given point of k-space anisotropy 
for both intraband and exchange interband terms. The 
intraband distributions of the two bands show different 
shapes and widths and have different range of values. 
The resulting matrix of coupling terms is obtained ex- 
clusively from the eigenfunctions of the superlattice and 
is asymmetric. 

Below, we present numerical results for the solution 
of the self-consistency equations which determines the 
values of the gap in the large FS and in the small FS. We 
discuss the behavior of Ai and A2, defined as average 
values of the gaps on the corresponding branches of the 
Fermi surface, and of the critical temperature T c . 

We discuss the case where the coupling term in the 
second band is smaller than in the first band (called here 
the pnictide-like case) . This case allows to emphasize the 
role of the exchange-like interband pairing and to use the 
BCS weak limit for interactions for contact interactions, 
neglecting retardation effects. In the extreme case when 
the intraband coupling in the second band is zero, the 
pairs in the second band are formed only thanks to in- 
terband exchange-like pairing. 

Therefore in this approach we consider the intraband 
coupling C1.1 for the band that does not show the Lifshitz 
transition, the intraband coupling 02,2 for the one that 
shows the Lifshitz transition and the inteband coupling 
Ci 2 for the exchange-like pairing. 

Fig. 4 shows the case of weak coupling in the second 
band where the intraband coupling parameters ratio is 
fixed at C2.2/C1.1 = 0.45. The critical temperature T c 
(panels a ), the ratio 2Ai/T c (panel b), and the ratio 
— 2/^.2 /T c (panel c) are plotted as function of the re- 
duced Lifshitz parameter (/x — E^lV^o. The different 
curves in each panel represent the cases of different in- 
terband pairing strengths (ci^/c^i =-0.68, -0.91, -1.04, 
-1.59, -2.73). In correspondence of the largest interband 
pairing strength examined here, C1.2 =-2.73, we have ob- 
tained critical temperatures as large as 50K, in the range 
1 < (/i — i?2,i)/wo < 2. The antiresonance in the first 
gap appears when the ratio 2Ai/T c reaches the minimum 
value, that can be much smaller than standard BCS value 
(2A/T C = 3.53). The antiresonance appears in the range 
— 1 < {n — E 2 ,l)/ l < j o < and moves from -1 to zero 
decreasing the interband pairing. Moreover the second 
superconducting gap values results to be nonzero even 
before the lower band edge E 2: l is reached. We notice 
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FIG. 5. (Color online). The effect of the variation of transver- 
sal dispersion on the shape resonance, the case of weak cou- 
pling in the second band (the so called iron pnictide case) for 
the case of strong interband repulsive interaction. The intra- 
band and interband coupling parameter ratio (02,2/01,1 = 0.45 
and ci, 2/01,1 = -3.6, with ci,i = 0.22) are fixed. The criti- 
cal temperature (panel a), the ratios 2Ai/T c (panel b) and 
— 2A-2/T c (panel c) as function of the reduced Lifshitz pa- 
rameter (fi ~ E2,l)/uo are plotted for £ = 110 (blue curve 
with open circles), 68 (red curve with open triangles) and 36 
meV (black curve with open squares). This last case is close 
to the case of a single quantum well. We also plot the DOS 
of the second subband in the inset of panel (a). The ratio 
2Ai/T c , probing the pairing channel in the first large FS, ex- 
hibits two minima due to the negative interference effect (or 
antiresonance) typical of shape resonances, at the band edge 
(M - E 2 ,l)/oj =0 and the the ETT (fx - E 2 , L )/uo=-l- The 
antiresonance effect in the large FS increases for decreasing 
dispersion in the second band. The value of the second su- 
perconducting gap is nonzero even before the lower band edge 
E2,l is reached. In the second band the maximum value of 
the ratio is always reached in the Bose regime and remains 
the same, — 2A2/T C ~ 6.7, for all dispersion values. 



that the maximum T c is reached in the BCS regime zone 
near the type II Lifshitz transition (ETTsd-2d)- 

Iron pnictides superconductors are certainly more 
three dimensional than cuprates (they have a quasi- 
2D electronic structure), but the dimensionality of elec- 
tronic structure in these compounds is also less clear, in 
fact even if some band structure works predict strong 
k 2 dispersions^ other works predict weak dispersions 
in AFe2As2 and ReOFeAsji^ determined mainly by Fe 
d-orbitals of FeAs layers. Furthermore many angle- 
resolved photoemission spectroscopy experimentsi22r-i2& 
reveal weak k 2 dependences for the bands electronic 
structure, in fact most of the band dispersion of 1111 
and 11 iron pnictides, like LaOFeAs and NaFeAs) 109 ' 110 
are in the range of 20 ~ 40 meV. In 122 undoped pnic- 
tides, like BaFe2As2, the k z dispersion is rather smalljiii 
although in BaFe2-xCo a; As2, with increasing Co concen- 
tration, the k z dispersion increases up to 120 meVJ^ 

Fig. 5 shows the effect of the variation of out of plane 
dispersion on the shape resonance, in the case of weak 
coupling in the second band (the so called iron pnictide 
case) for strong interband pairing. The DOS of the sec- 
ond subband for different dispersions (£ = 110, 68 and 
36 meV), as a function of the reduced Lifshitz param- 
eter (fi — £^2,l)/wo, is shown as an inset in panel (a). 
The intraband and interband coupling parameter ratio 
( c 2,2/ci,i = 0.45 and ci, 2/01,1 = -3.6, with ci,i = 0.22) 
are fixed. The critical temperature, panel (b), the ra- 
tios 2Ai/T c , panel (c), and —2A 2 /T C , panel (d), are 
plotted as a function of the reduced Lifshitz parameter 
(/i — E2,l)/ uj o- The minima of the ratio 2Ai/T c , prob- 
ing the pairing channel in the first large Fermi surface 
show two minima due to the negative interference effect 
or antiresonance typical of shape resonances at the band 
edge (n - E 2 x)/lj =0, and the ETT (/x - E 2tL )/u =-l. 
The antiresonance effect in the large FS increases for de- 
creasing dispersion in the second band. The value of the 
second superconducting gap is nonzero even before the 
lower band edge E 2 .l is reached. In the second band, 
the maximum value of the ratio is always reached in the 
Bose regime and remains the same, —2A 2 /T C ~ 6.7, for 
all dispersion values. 

Using STM spectroscopy it is now possible to measure 
the gaps and the critical temperature in the same set 
of experiments. Therefore the BCS gap ratios 2Ai/T c 
and -2A 2 /T C as functions of the ratio T c /T™ ax can be 
measured directly for different samples of different gate 
voltages tuning the chemical potential. The present the- 
ory is able to predict these curves and we show that these 
curves provide a direct measure of the relevance of inter- 
band coupling versus Cooper pairing. In fig. 6 we show 
the case strong interband repulsive interaction and weak 
coupling for Cooper pairing. We keep constant the ra- 
tio C2, 2/01,1 = 0.45, ci.i = 0.22 and £/ojo = 1 and we 
present the expected behavior of these curves for strong 
interband interaction. The BCS Gap ratio parameters 
are plotted for three cases of different interband coupling 
parameters. When the interband repulsive term is domi- 
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FIG. 6. (Color online), (a) The ratio 2Ai/T c and -2A 2 /T C as 
function of the ratio T c /T™ ax , in the case of weak coupling in 
the second band, with fixed 02,2/01,1 = 0.45, ci,i = 0.22 and 
interlayer dispersion (£/ojo = 1), for three cases of different 
interband coupling parameters (01,2/01,1 = -2.27, -3-63, and 
-5.45) ratio, (b) The ratios 2Ai/T c and — 2A 2 /T C as functions 
of T c in the case of weak coupling in the second band with 
fixed 02,2/01,1 = 0.45, ci,i = 0.22 and interlayer dispersion 
(£/wo = 1), and in the case of weak coupling in the second 
band with fixed 01,2/01,1 = 0.45, ci,i = 0.22 and interlayer 
dispersion £ = 110 (blue curve with open circles), 68 (red 
curve with open triangles) and 36 meV (black curve with open 
squares) . 

nant, the BCS gap ratio is the same in both bands at the 
maximum critical temperature. The difference between 
the gap ratios diverges decreasing the critical tempera- 
ture. This behavior is very similar to case of cuprates in 
fact in the underdoped region where the critical tempera- 
ture goes to zero and while a first BCS-like gap decreases 
with the decreasing critical temperature, the second gap 
(called pseudogap) increases by decreasing the critical 
temperature. 

We have investigated the effect of the band dispersion 
for a case on moderate relevance of the interband pairing 
such that there is a characteristic critical temperature 
where the two gaps are the same and near the standard 
BCS gap parameter 3.53. We can see that going from 
the case where the band dispersion is low, i.e., we are 



close to the regime for a single isolated slab of Blatt to 
the regime where the dispersion is as large as twice the 
energy cutoff of the interaction there minor variations. 
This result shows that the shape resonance due to quan- 
tum confinement is a robust feature also for appreciable 
electron hopping between the quantum wells. 



VI. CONCLUSIONS 

The present work provides an interpretation for the 
properties of quantum size effects in superconducting 
multilayers and indicate a possible roadmap for the dis- 
covery of novel HTS like graphene bilayers and graphane 
superlattices that should share similarities with the 
known multigap HTS families where the shape resonance 
is driven by the interband pairing mechanisms. 

Here, we have investigated the effect of electron hop- 
ping between quantum wells that can be changed by 
changing the thickness of spacer layers, like in 1111 pnic- 
tides where the rare earth ionic radius in the spacer layer 
can be changed, or by changing the type of material form- 
ing the spacers like going from 1111 to 122 pnictides. To 
take into account this effect on the shape resonance, we 
changed the dispersion £, while keeping it of the order of 
the pairing interaction energy cut-off u>q. We have fixed 
in our calculation the intraband Cooper pairing much 
weaker than the exchange-like interband repulsive inter- 
action. This is clearly the case of iron pnictides that 
exhibit critical temperatures of the order of 50K mostly 
driven by interband pairing. Namely, we have examined 
the case of intraband and interband coupling parameter 
ratio C2,2/ci,i = 0.45 with variable interband repulsive 
interactions and variable single electron dispersion in the 
transversal direction. 

The ratios 2Ai/T c and -2A 2 /T C predicted for high-T c 
iron-based layered compounds and cuprates as functions 
of T c show that for the highest T c samples the gaps be- 
come equivalent while they arc quite different in the low 
T c samples, that appears to be the scenario for cuprates 
and for pnictides. The antiresonance typical of shape res- 
onance occurs where the Lifshitz energy parameter is zero 
(i.e appearing of a new Fermi surface spot) for weak in- 
terband exchange pairing and shift to -1 for a very strong 
interband exchange pairing. The relevant result is that 
the maximum of the critical temperature appears near 
the Lifshitz transition of the type opening a neck in cor- 
rugated cylindrical surface therefore it moves with the 
variation of the electronic hopping between the layers i.e., 
the energy dispersion in the normal direction respect to 
the superconducting layers. 
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